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Abstract —The growing amount of fluctuating renewable in¬ 
feeds and market liberalization increases uncertainty in power 
system operation. To capture the influence of fluctuations in 
operational planning, we model the forecast errors of the un¬ 
certain in-feeds as random variables and formulate a security 
constrained optimal power flow using chance constraints. The 
chance constraints limit the probability of violations of technical 
constraints, such as generation and transmission limits, but 
require a tractable reformulation. In this paper, we discuss 
different analytical reformulations of the chance constraints, 
based on a given set of assumptions concerning the forecast 
error distributions. In particular, we discuss reformulations that 
do not assume a normal distribution, and admit an analytical 
reformulation given only a mean vector and covariance matrix. 
We illustrate our method with a case study of the IEEE 118 
bus system, based on real data from the European system. The 
different reformulations are compared in terms of both achieved 
empirical violation probability and operational cost, which allows 
us to provide a suggestion for the most appropriate reformulation 
in an optimal power flow setting. For a large number of 
uncertainty sources, it is observed that the distributions of the 
line flows and generator outputs can be close to normal, even 
though the power injections are not normally distributed. 

Index Terms —Renewable Integration, Chance Constrained 
Optimal Power Flow, N-I security 

I. Introduction 

A fundamental tool in power system analysis is the optimal 
power flow (OPF) IT]. Several tasks central to power system 
operation, such as unit commitment, reserve procurement, 
market clearing and security assessment rely on the solution of 
an OPF. The main goal of the OPF is to minimize operational 
cost, while ensuring secure operation that respects technical 
limits of the power system. In current operational schemes, 
the system is considered secure if it remains within the 
operational limits during normal operation and during outage 
of any single component. This principle is referred to as the 
A^ — 1 criterion, and is reflected in the OPF through additional 
constraints, leading to a security constrained optimal power 
flow (SCOPF). While the N — 1 criterion secures the system 
against individual outages,/orecavf uncertainty is another kind 
of disturbance affecting the system. Forecast uncertainty arises 
from unforeseen fluctuations in the power injections, such as 
inaccurate predictions of load or renewable in-feeds, as well 
as from short-term electricity trading. While load profiles are 
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relatively predictable, higher shares of electricity production 
from renewable sources and liberalization of energy markets 
(particularly in Europe) have increased the forecast uncertainty 
by orders of magnitude El. In current operational planning, 
uncertainty is usually ignored and uncertain quantities are 
typically replaced by a forecast value. While this approach has 
provided good solutions in the past, the increased levels of un¬ 
certainty lead to frequent N-1 violations in real-time operation. 
To mitigate these problems, it is proposed to explicitly account 
for uncertainty during operational planning, in particular while 
solving the OPF. 

There are different approaches to account for uncertainty 
within the OPF. Robust and worst-case methods, e.g. E), en¬ 
sure secure operations for all possible forecast errors, but often 
provide very conservative and thus costly solutions. Stochastic 
programming methods give the operator more freedom to 
trade-off cost and security. One example is two stage stochastic 
programming for unit commitment and reserve scheduling, 
e.g. a, which minimizes the expected cost of operations 
based on a set of scenarios. Another example is chance con¬ 
strained programming, which explicitly limits the probability 
of constraint violations a, a, a. Since the main goal in 
short-term operational planning is to ensure secure operations, 
we consider the latter method, and formulate the OPF as a 
probabilistic SCOPF (pSCOPF) with chance constraints. The 
acceptable violation probability, which is treated as a design 
parameter in the optimization problem, allows the operator to 
choose an appropriate trade-off between cost and security of 
operations. 

Although the pSCOPF allow us to account for uncertainty 
in a comprehensive way, it is generally hard to reformu¬ 
late chance constraints as tractable constraints. Two main 
approaches for reformulation have been applied to the OPF 
problem, based either on sampling or analytical reformulation. 
In |5l, the SCOPF is formulated as a joint chance constrained 
problem (limiting the probability that any of the constraints are 
violated), which is reformulated using the scenario approach 
based on 0. The formulation was extended to include market 
clearing with co-optimization of energy and reserves in Q, 
where a different sampling based reformulation based on 
HI was used. Both sampling based reformulations require 
no knowledge about the underlying distribution, except for 
availability of a given number of samples (which increases 
with the problem size). 

In contrast, the SCOPF formulated in Q uses separate 
chance constraints (limiting the probability for each constraint 
separately) to formulate the pSCOPF. Assuming that the 
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random variables follow a Gaussian distribution, an exact 
analytical reformulation is obtained. The same type of Gaus¬ 
sian reformulation is performed for an OFF without security 
constraints in la. 

While the assumption of a Gaussian distribution limits the 
applicability of the analytical reformulation from Q, (a, the 
analytical reformulation has some attractive properties. First, 
it is scalable to a large number of random variables, as the 
number of random variables does not influence the problem 
size or complexity of the OFF itself. Second, the solution 
is more transparent than a sample based solution since it 
is possible to trace the influence of each random variable 
through the analytical relations. Finally, the solution based on 
the analytical reformulation is deterministic, i.e., the OFF will 
always And the same optimal solution with the same optimal 
cost. While this might seem trivial, the OFF solution based 
on the scenario approach is actually random, since it depends 
on the choice of the samples. The same problem might thus 
lead to different solutions with different costs, depending on 
which samples were chosen. 

This paper investigates how the good qualities of the an¬ 
alytical reformulation can be preserved, while moving away 
from the limiting assumption of a Gaussian distribution. Using 
optimal probability inequalities, we obtain distributionally 
robust reformulations of the chance constraints. This approach 
is well-known in operations research and control theory and 
has been investigated in, e.g., ifTTI . ifT^ . lfT3l . In lfT4l . the 
application of some distributionally robust reformulations to 
the optimal power flow problem were also discussed. Here, 
we introduce reformulations based on assumptions like uni¬ 
modality and symmetry of the forecast errors, and explain why 
those are relevant in the optimal power flow context. We aim to 
provide recommendations for the most suitable reformulations, 
depending on the sources of uncertainty (e.g., RES fluctua¬ 
tions, load variations, short term trading) and the time frame 
(e.g., day-ahead planning, real-time operation). To compare the 
different analytical reformulations, we introduce the concept of 
an uncertainty margin. The uncertainty margin has a physical 
interpretation as a security margin against forecast errors, and 
represents a reduction of available transmission and generation 
capacity. A larger uncertainty margin thus increases security, 
but also the operational cost. The empirical performance of 
the proposed reformulations is assessed through a case study 
for the IEEE 118 bus system, with uncertainties represented 
through historical forecast errors from the Austrian Fower 
Grid. We investigate which reformulation is the most appro¬ 
priate for a chance constrained SCOFF problem, considering 
empirical violation probability, nominal operational cost and 
the accuracy of the distributional assumptions. 

The remainder of this paper is organized as follows: Section 
II and III present the uncertainty modeling and the formulation 
of the chance constrained SCOFF. Section IV discusses differ¬ 
ent analytical reformulations for these chance constraints, and 
Section V applies them to the power flow equations. Section V 
demonstrates the proposed formulation in a case study for the 
IEEE 118 bus system. Section VI summarizes and concludes 
the paper. 


II. Modeling forecast uncertainty 

Forecast errors arise from uncontrolled power in-feeds that 
deviate from their forecasted values, caused by, e.g., fluctua¬ 
tions in load or renewable production. The characteristics of 
forecast errors Sr differ between systems, depending on their 
generation mix, load characteristics and market structure, and 
also depend on parameters such as the time of the day, or the 
forecast horizon. For example, the minute-to-minute variation 
in wind in-feeds in the central European system follows a 
Student t-distribution CSl, while the day-ahead distribution 
of wind forecast depends on the forecasted wind power, 
and is typically non-symmetric. The uncertain in-feeds can 
be modeled as random variables with continuous probability 
distributions. We define the vector of n uncertain in-feeds as 

Pr = Pr + Sr. (1) 

Here, the uncertain in-feed Pr is the sum of the forecasted 
value Pr and a random deviation Sr. Since the sources of 
uncertainty differ both within and between systems, the full 
distribution of Sr is generally not known. However, we will 
assume some partial information about the distribution. In 
particular, we assume that the mean fXR S K" and covariance 
PiR S of the forecast errors exist, and can be estimated 

either based on historical data or through forecasting methods. 
We allow for non-zero mean, since forecasts are not necessar¬ 
ily based on the expectation of Pr, but rather on the most 
probable realization (which are not the same, e.g., for skewed 
distributions). 

HI. Optimal Fower Flow Formulation 

We now introduce the mathematical formulation of the 
pSCOFF for a system with n buses, based on the formulation 
in Q. The sets Q, T) and TZ represent the conventional 
generators, the fixed loads and the uncertain in-feeds (con¬ 
sisting of, e.g., in-feeds from wind and solar power plants), 
respectively. To simplify notation, we assume that there is 
one generator, load and generator connected at each bus, such 
that \Q\ = \V\ = \R.\ = n. This assumption is however not 
necessary for the method itself. The set of transmission lines 
is denoted by C, and there are |£| = ul lines in the system. 
The contingencies considered for the N-1 security criterion 
include outage of any line or generator, in total nc = ul +n 
outages. 

A. Generator modeling 

The nominal generation output of the generators, Pq € K." 
are the optimization variables of the problem. In addition to 
keeping the system balanced in nominal operation conditions, 
any power deviation arising from either forecast errors or 
generation outages must be balanced by the generators. The 
contribution of balancing energy from each generator can be 
chosen in different ways. Here, we assume that each generator 
contributes according to its maximum nominal output, similar 
to |[T6l . When all generators operate, the balancing contribu¬ 
tion of each generator g is given by 

pmax 

gi — ^’3 I'T'v 

“(s) pmax ’ t 
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where the superscript 0 refers to normal operating condition 
{i = 0), or situation with line outages {i G C). During 
the outage of generator i, the compensation vector of the 
generators is given by 

Vpv, = 0 . (3) 

The vectors d® G R” thus describe the compensation of any 
power mismatch in the system, for any outage situation i. We 
note that by definition of d, the system remains balanced after 
any fluctuations or generator outages. 


dU = 


jDmax 

^G,g 


E n 

i=i.. 


pmax 


B. Line flow modeling 

Similar to the setup in 0, the line flows are expressed as 
linear functions of the active power injections in both normal 
and outage conditions. 


Pi = A^Pl^ for alH = 0,...,nc . 


(4) 


Here, A® G R®®r,xn describes the relation between the active 
power injections P^^^ G K®® and the line flows PJ after outage 
i, with i = 0 being the normal operation condition. A® is given 

by 


A^ = B} 



0 

0 


(5) 


depend on the realization of the random variable Sr, and are 
formulated as single chance constraints. The chance constraint 
ensures that probability of a constraint violation (e.g., a line 
flow exceeding the limit) remains smaller than e. We will refer 
to e as the violation probability and to 1 — e as the security 
level. The value of e is an input parameter to the optimization. 

IV. Chance constraint reformulation 

To obtain a tractable optimization problem, the chance 
constraints must admit a deterministic and tractable 

reformulation. These constraints (l9ll- (fT2l i are all univariate or 
single chance constraints of the general form 

P[a(PG) + biPG)SR < c] > 1 - e . (13) 

where a^Pc) G R and b{PG) G are affine functions of 

the decision variables Pq and c is a constant. The term a{PG) 
represents the nominal generation output or the nominal line 
flows (without forecast errors) and c represents the generation 
or line flow limit. The vector 6 (Pg) expresses the influence of 
the forecast errors 6 r on the respective constraint. Regardless 
of the exact expressions for b{PG), and for any dimension or 
distribution of the random vector 5 r, the left hand side of the 
constraint is a scalar random variable 5 = a{PG) + b{PG)5R 
with mean h^Pq) and variance ct^Pg) given by 


where Pj G M"r.xn jjj^g sus^eptance matrix and 

^bus ^ the bus susceptance matrix (without the 

last column and row) after outage i 0. The power injections 
are given by 

^inj = -Pg + d"(PG(i) - Is^r) + Pr + Sr - Pr . (6) 

Here, Pr G R®® is the vector of loads, and the power mismatch 
due to generation outage PG(i) is non-zero only for i G Q. 
The vector 1^ G R^®^®® is a vector of ones, such that 
represents the sum of the forecast errors. 

C. Chance constrained optimal power flow 

Using the modeling assumptions presented above, we can 


formulate the pSCOPF as 

min (F Pg (7) 

Pc 

subject to 

llxn(PG + Pr — Pl) = 0 (8) 

nPGig) + (PG(i) - USr) < P5®(“f ] > 1 - £ , (9) 

nPGig) + rf(s)(PG(.) - IsSr) > PS®(®;5] > 1 - £ , (10) 

nA\i,.)PLj<Pml>^-e, (11) 

> -Pi^f] > 1 - ^, (12) 

for g = 1, I = 1, ...,nL, i = 1, ...,nc • 


The objective 0 is to minimize generation cost, with c 
representing the bids of the generators. Constraint ® ensures 
power balance in the system. The constraints (l9li- (fT2li are the 
generation and transmission constraints, with p™™ and 
being the minimum and maximum generation levels and 
being the transmission capacity of the lines. Those constraints 


fJ-iPG) = ^{Pg) + b{PG)p-R , O'(Pg) =11 ^(Pg)^^^ II2 • 

What is of interest when reformulating the constraint (fT3l l is 
not the distribution of the forecast uncertainty 6r, but the 
distribution of S, which represents the variations in line flows 
or generation outputs. Depending on the system, S might 
follow different distributions. We will now present different 
distributional assumptions for S which are relevant in the 
context of the SCOPE The applicability of each assumption 
depends mainly on the source of uncertainty (e.g., load, renew¬ 
ables or short-term trading), the time frame of the forecast 
(e.g., day-ahead planning or close to real-time operation) 
and the availability of data (e.g., historical forecast errors or 
probabilistic forecasts). 

1) Normal distribution ($).' The normal distribution is a 
good distribution model in two different cases. First, when 
5 r follows a multivariate normal distribution (which might 
be the case, e.g., for load uncertainty), which means 5 will 
be normally distributed as well. Second, when the number 
of uncertainty sources is large and not highly correlated, 
arguments similar to the central limit theorem (e.g., CD) 
imply that the distribution of <5 (which is a weighted sum of 
5 r) is expected to be close to a normal distribution. 

2) Student’s t-distribution (t): When the forecast fluctua¬ 
tions are heavy tailed (e.g., as for the minute-to-minute vari¬ 
ability in the European grid ifTSl l. the Student’s f-distribution 
can be a more appropriate representation. Particularly when 
considering small violation probabilities (e < 0.03), Student’s 
t distribution provides additional robustness compared to the 
normal distribution. 

In many cases, only limited knowledge about the distri¬ 
bution of 6 is available. It might therefore be desirable to 
only assume some general properties of the distribution of S, 
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rather than a specific distribution. This leads to the following 
distributionally robust reformulations, that are valid for all 
probability distributions that share the general properties: 

3) Symmetric, unimodal distributions (S): If the distribu¬ 
tion is likely to be close to normal, but we do not know how 
close, we can resort to the general assumption of unimodal, 
symmetric distribution with known mean and covariance. 

4) Unimodal distributions (U): In systems where the fore¬ 
cast uncertainty is related mainly to load, wind and PV 
production, the distribution of Sji is likely to be unimodal, 
with fluctuations centered around the forecasted value. Under 
such conditions, it is highly probable that the distribution of 
6 is also unimodal. 

5) Known mean and covariance (C): In systems where 
intra-day electricity trading is not controlled by the trans¬ 
mission system operator, for example in Europe, intra-day 
transactions introduce uncertainty in the power injections from 
conventional power plants. The transactions might follow 
almost any probability distribution, and can even be discrete. 
In this case, we reformulate the chance constraint based only 
on a known (and finite) mean and covariance. 

For all distributional assumptions 1) - 5), the chance con¬ 
straint O can be reformulated to the following analytic 
expression 

a{PG)<c-b{PG)tiR-r\^-e) II KPg)^T II 2 ■ (14) 

Analyzing (O, we see that the left part, a{PG) < c, represents 
the “nominal” constraint, i.e., the constraint we would obtain 
if we neglect the forecast uncertainty. The second and third 
term represents a reduction of the nominally available capacity 
c, which is necessary to secure the system against forecast 
deviations. This reduction can thus be interpreted as a security 
margin against uncertainty, i.e., an uncertainty margin. Notice 
that the larger /“^(l — e), the larger the uncertainty margin. 

Depending on which assumption 1) - 5) is deemed appro¬ 
priate, we define /~^(1 — e) according to either an inverse 
cumulative distribution function (for known distributions 1), 
2)) or a probability inequality (when only partial information 
is available 3) - 5)). The exact expressions for /“^(l — e) are 
shown in Table |I] and their derivations as well as the derivation 
of (fT4l) are given in the Appendix. We note that for 1) and 2), 
the reformulation is tight (the chance constraint holds with 
equality). The distributionally robust reformulations 3)-5) are 
typically not tight, and will usually lead to empirical violation 
probabilities lower than e. 

Since the reformulations 1) - 5) differ only in the definition 
of /~^(1 — e), we can compare them by comparing the value 
of /“^(l — e) for different e. In Fig. [T] /“^(l — e) is plotted 
against the security level 1 — e. We observe that all /“^(l — e) 
increase as e decreases, indicating that a larger uncertainty 
margin is necessary to achieve a lower violation probability. 
With more information, tighter probabilistic bounds can be 
defined and thus a lower value of /“^(l — e) is necessary 
to ensure the desired security level (i.e., ^ > fg^). 

The lowest values are obtained when we assume knowledge of 
the actual distribution, i.e., for the normal and the Student’s t 
distribution. Note that all reformulations assuming symmetry 
have fs{0.5) = /<i>(0.5) = /t(0.5) = 0. Finally, Student’s 


TABLE I 

Expressions for /“^(l — e). 

Cumulative distribution function of the standard normal distribution. 
Cumulative distribution function of the Student t distribution with 
zero mean, ly degrees of freedom and scale parameter ctt — (^ ~ 2)/r/. 


1) Normal 

1 

1 

II 

1 

1 101 

2) Student’s t 


3) Symmetric, unimodal 

/s(l - e) = • 

' 

VT for 0 < £ < 1 

\/3(l — 2s) for g < £ < 5 

.0 for 4 < £ < 1 

4) Unimodal 

/c7^(l-e) = • 

V ^ - 1 for 0 < £ < 4 
, V l+3e for 1 < e < 1 

5) Mean, covariance 

= forO<e<l 




Security level 1-e 


-fJ Norma! 

<t> 

t with v=4 

-Chebyshev 

Unimodal 

-Symm-fUni 


Fig. 1. Values of /”^(1 — e) for the normal distribution, the Student t 
distribution with 4 degrees of freedom, the Chebyshev inequality, unimodal 
distributions and symmetric, unimodal distributions. The left part shows all 
security levels, while the right part is a zoom in on high security levels. 


t distribution has a more pronounced peak and heavier tails 
than the normal distribution. This is reflected in that for lower 
security levels, while at high security 

levels, /”^(1 - e) > /$ ^(1 - e). 

V. Reformulated constraints 

With the reformulation presented above, we can reformulate 
the chance constraints d^- lfTSI) as 

PGi,)+dl^PGi^)<P^C,)- (15) 

dl^lsf^R-r\l-e) lb, 

PGi,)+dl)PGii)>P^l-)- (16) 

d(,)l5Mfl + r'(l-e) lb, 

(Pg + d^PGH) +Pr- Pd) < P^d)- (17) 

A^^i,)[Pg + ddPGU) +Pr- Pd) > -P^df- (18) 

+ /-i(l-e) II lb, 

for 5 = 1, ...,n, I = 1, ...,n, i = 1, ...,nc , 

where / g is the identity matrix. When comparing ( fTSl ) 

- (O to (114b . we recognize the same structure. The first part 
represent the constraint we would obtain by neglecting the 
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forecast uncertainty, where as the second and third term on 
the right hand side represent the uncertainty margin. 

Since a higher uncertainty margin leads to a reduction in 
the available transmission and generation capacity, a higher 
uncertainty margin will not only reduce the probability of 
violation, but also increase the nominal cost of operation (i.e., 
the cost of the pSCOPF). The acceptable violation probability 
£ and the distributional assumption (which defines the function 
f~^{l — e)) should therefore be chosen carefully to obtain a 
good trade-off between security against forecast errors and 
cost of operation. 

Note that the reformulated chance constraints dnil - (HI 
are linear, since the uncertainty margin is not dependent on 
any decision variables and can be pre-computed. The pSCOPF 
problem (|7]i, (IHl, (flST l - (fTsT l is thus a linear program with the 
same computational complexity as a traditional DC SCOPF. 

VI. Case Studies 

The purpose of this case study is to demonstrate the chance- 
constrained SCOPF, and investigate which distributional as¬ 
sumptions are most appropriate for power systems operation. 
We base our study on the IEEE 118-bus system M, with a 
few modifications as follows. The generation cost is assumed 
to be linear, and is based on the linear cost coefficients 
of the data provided with Matpower 4.1 QS). Although the 
formulation could be extended to include unit commitment, it 
is not considered here. Therefore, the minimum generation 
output of the conventional generators is set to zero. The 
forecast uncertainty Sr is modeled based on historical data 
for 1 year from the Austrian Power Grid (APG). We define 
the forecast error as the difference between the the so-called 
DACE (Day-Ahead Congestion Eorecast) and the snapshot 
(the real-time power injections) for all hours and buses with 
available data (8492 data points for 28 buses). Since the system 
is constantly evolving and might exhibit seasonal patterns, we 
assume that the power system operator only uses data from the 
past three months. We use two three-month periods to define 
the forecast uncertainty for this case study, such that we obtain 
2207 data samples for a total of 54 buses. 

The historical data was assigned to different load buses 
throughout the system, and modified such that the standard 
deviation corresponds to 20 % of the forecasted load. The 
mean and covariance E/j used in the pSCOPE were 
calculated based on this modified data (i.e, assuming perfect 
knowledge of E^;). Eig.|2] shows the forecast errors from 
some representative nodes, including the histograms and pair¬ 
wise scatter plots of the forecast errors. By inspection, it is 
clear that the forecast errors are not normally distributed. 

In the following, we assess how the different distributional 
assumptions impact the solution of the pSCOPE. We solve the 
pSCOPE for all five reformulations assuming an acceptable 
violation probability of e = 0.1. The results are compared 
with each other and to the solution of the corresponding 
deterministic SCOPE. To assess the quality of the solution, we 
compare the number of empirical constraint violations (based 
on the historical samples) and the relative cost of the solutions. 
Eurther, we run statistical tests to check if the data is normally 
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Fig. 2. Forecast errors for 4 selected nodes of case study. The diagonal plots 
show the histograms of the forecast eiTors (x-axis: deviation in MW, y-axis: 
number of occurences), while the off-diagonal plots show the scatter plots 
between two corresponding forecast errors (x- and y-axis: deviation in MW). 
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Fig. 3. Results derived from the different SCOPF solutions. From left to right: 
0) the deterministic SCOPF, and the pSCOPF based on the assumption of 1) 
a normal distribution, 2) a Student t distribution, 3) a symmetric, unimodal 
distribution, 4) a unimodal distribution and 5) known mean and covariance. 
From top to bottom, the figure shows a) empirical violation probabilities for 
the all constraints, b) nominal dispatch cost, and c) value of the function 
/“^(l — e) for £ = 0.1. 


or unimodally distributed, and investigate the accuracy of the 
estimated uncertainty margins. 

1) Number of empirical violations: The empirical violation 
probabilities e are evaluated for all constraints based on the 
2207 data samples. The results are shown in Eig. [3a), starting 
with 0) the deterministic solution, then the solution based 
on 1) a normal and 2) a Student t distribution, and then 
the distributionally robust solutions 3)-5). Erom left to right, 
we thus assume lessened knowledge about the distribution. 
The violation probabilities of the non-active constraints in-a 
and active ia constraints are plotted in yellow and orange, 
respectively. The average violation probability iavg of the 
active constraints is plotted in black. 

As seen on the left, the empirical violation probability of 
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the deterministic solution is very high, with some constraints 
violating 80% of the cases. This highlights the need for prob¬ 
abilistic methods to avoid frequent violations of operational 
limits. The pSCOPF solutions have much lower violation 
probabilities. The active constraints (where the distance from 
the nominal flow to the flow limit is given by the uncertainty 
margin) have a higher empirical violation probability than the 
non-active constraints (which have some additional margin). 
The solution based on a normal distribution violates the 
accepted violation probability e < 0.1 for some constraints, 
but the violation is small, e < 0.11, and the average vi¬ 
olation probability is acceptable iavg < 0.1. The solution 
based on a Student t distribution, which assumes a more 
peaked distribution, has larger empirical violation probabilities 
e > 0.15. The distributionally robust solutions oversatisfy 
the accepted violation probability e, with iavg = 0.07 for 
the symmetric, unimodal solution, 0.03 for the unimodal 
solution and 0.0025 for the solution based only on mean and 
covariance. Although the chance constraints are satisfied, it 
does not necessarily imply that the underlying assumption 
(e.g., symmetry and unimodality) is accurate. Since the re¬ 
formulations are distributionally robust, we might get a low 
empirical violation probability, even if we assumed the wrong 
family of distributions. 

2 ) Operational cost: Fig [3] b) shows the generation cost 
obtained with the pSCOPF solutions, normalized by the cost of 
the deterministic problem (shown to the left). All probabilistic 
solutions have higher cost than the deterministic solution, 
showing that the consideration of uncertainties increase the 
nominal cost of operation. The reformulations which assume 
more knowledge about the distribution 1), 2) lead to lower 
cost than the more general reformulations 3)-5). The most 
expensive solution is obtained for reformulation 5), which only 
assumes knowledge of mean and covariance. 

The cost differences are explained by the different values 
— which defines the uncertainty margin and thus the 
constraint tightening. The value of f~^{l—e) is plotted in Fig. 
[3c). Comparing Fig.|3]a), b) and c), we observe how a larger 
/“^(l — e) leads to an increase in nominal operation cost, but 
at the same time reduces the empirical violation probability. 
This highlights two important aspects of the pSCOPF. First, 
we need to define e such that it reflects a reasonable trade¬ 
off between cost and security. Second, we want to achieve 
an empirical violation probability e as close as possible to 
the accepted violation level e. A reformulation with too many 
violations (i >> e) leads to unsecure operations, but at the 
same time a too conservative solution (e << e) will lead to 
unnecessary high cost, and possibly infeasibility if we want to 
ensure a low violation probability (e.g., e ss 0.01). 

3) Testing the distributional assumptions and the accuracy 
of estimated uncertainty margin: Since the transmission and 
generation constraints are enforced as separate chance con¬ 
straints with deviations S (defined as a weighted sum of the 
random variables Sr), each constraint has a univariate distri¬ 
bution related to it. To assess whether or not our assumptions 
about those distributions are correct, we run statistical tests. 
In particular, we use the Shapiro-Wilk test li20l to test if 
the distribution is normal, and Hartigans dip test ED to test 


80 r 
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I Shapiro-Wilk Test for Normality 
I Hartigans Dip Test for Unimodality 
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P-values of the Shapiro-Wilk and Hartican’s Dip Test for all constraints 


Fig. 4. P-values obtained from Hartigans dip test for unimodality and the 
Shapiro-Wilk test of normality. The histogram show the percentage of all 
constraints in each 0.05 interval. A high p-value indicate a high probability 
that the distribution is unimodal or normal, respectively. Based on the results, 
it seems highly probable that the distributions are unimodal (p-values close 
to 1), while it is unlikely that the data comes from a normal distribution 
(p-values close to 0). 


unimodality, using the implementation in R li22ll . The test 
output is a p-value between 0 and 1, which indicates how 
probable it is that the data comes from a normal or a unimodal 
distribution, respectively. Typically, the hypothesis (normality 
or unimodality) is accepted for p-values above p > 0.95, and 
rejected for p-values p < 0.05. In between, we can neither 
reject nor confirm the hypothesis. 

In Fig. m the p-values from both tests are plotted as a 
histogram. The bars show the percentage of constraints with 
p-values in the indicated p-value interval. We observe that the 
p-values from the Shapiro-Wilk test are below the threshold 
p < 0.05 for most constraints, while the p-values from 
Hartigans Dip Test are above p > 0.95 for the majority of the 
constraints. We thus conclude that unimodality is a reasonable 
assumption, while the original data is probably not normally 
distributed. 

Although the statistical test rejects normality, the normal 
distribution might still be a good assumptions for the parts 
of the distribution which we are interested in. In Fig. [5] 
the empirical distribution of the line flow deviations for one 
active transmission constraint is shown. This constraint had 
the lowest p-value among the active constraints in the Shapiro- 
Wilk test. Fig. [5] also show the empirical uncertainty margin 
(corresponding to an empirical violation probability e = 0.1), 
as well as the uncertainty margins obtained with the five 
distributional assumptions. We observe that the uncertainty 
margins based on the normal distribution (plotted in red) 
match very closely to the empirical margins (plotted in green). 
The Student t distribution underestimates the margin and the 
distributionally robust reformulations lead to too high margins. 

4) Discussion: Based on the above results, we conclude 
that the reformulation based on a normal distribution can 
provide a good trade-off between cost and security, partic¬ 
ularly in systems with a large number of uncertainty sources. 
Although assuming a normal distribution does not guarantee 
empirical violation probabilities e < e, the assumption might 
be useful when e can be interpreted as a guideline rather 
than a hard limit. Statistical tests and assessments as in Fig. 
13 can be used to assess whether the normal distribution is 
a reasonable approximation. If the system operator wants a 
higher confidence in enforcing the actual violation level e 
and is willing to tolerate a larger increase in operational cost. 
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Fig. 5. Histogram of post-contingency line flow deviations. The uncertainty 
margins are computed empirically (green), and for 1) a normal distribution 
(red), 2) a Student t distribution (yellow), 3) a symmatric, unimodal distribu¬ 
tion (magenta), 4) a unimodal distribution (light blue) and 5) a distribution 
where only the mean and covariance are known (dark blue). 


assuming a unimodal distribution would be reasonable. 

VII. Summary and Conclusion 

This paper discusses different analytic reformulations for 
chance constraints and their applicability in the pSCOPF 
context. The chance constraints are reformulated either by 
assuming a known probability distribution (such as normal 
or Student t distribution) or by using distributionally robust 
reformulations assuming general properties of the distribution 
(i.e., known mean and variance, symmetry, unimodality). 

The reformulated chance constraints all have a similar 
form, and are easily comparable. They are similar to the 
nominal constraints of the deterministic problem, except for 
the uncertainty margin (a security margin against forecast 
deviations), which represents a reduction of the transmission 
or generation capacity. With a larger uncertainty margin, the 
probability of violations decreases, but the nominal operational 
cost increases. Therefore, it is desirable to find a reformulation 
which leads to an uncertainty margin which is sufficiently 
large, yet as small as possible. 

In the case study based on the IEEE 118 bus system and 
forecast errors from Austria, the trade-off between security 
and cost is highlighted. Although the choice of reformulation 
differs between systems with different uncertainty charac¬ 
teristics, we show that the normal distribution might be a 
good approximation in cases where the acceptable violation 
probability can be interpreted as a guideline, rather than a 
hard constraint. If the transmission system operator wants to 
enforce the violation probability as a strict limit, choosing a 
more conservative, distributionally robust reformulations based 
on, e.g., unimodality will provide more confidence. 

In general, we believe that the pSCOPE with analytically 
reformulated chance constraints provides a transparent and 
scalable approach to assess the effect of uncertainty in power 
system operational planning. Euture work will investigate 
how the approach can be extended towards corrective control 
actions for uncertainties (e.g., HVDC and PSTs), to further 
reduce the cost of handling uncertainty. Eurther, we plan to 
investigate extensions towards AC power flow. 
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Appendix 


A. Exact tractable reformulations of single chance constraints 
We discuss here how the chance constraint O can be 
reformulated as the deterministic constraint (fT4l) . Using the 
equality 5 = a{Po) + b{Po)5}i it is clear that constraint (flTl) 
is equivalent to 

P[5 < c] > 1 - e (19) 


where the constraint (fT^ should be satisfied for all distribu¬ 
tions P consistent with the distributional assumptions made on 
5. The constraint (fT^ can equivalently be represented as 


S- tiiPg) c- fiiPc) 
(t(Pg) cr(-PG) 


> 1 - e 


( 20 ) 


where it can be remarked that the scaled random variable 
5n ■= {5 — ij,{Pg)) /a{PG) has zero mean and unit variance 
by construction. 

In order to unify the analysis for all distributional assump¬ 
tions made in Section HVl we consider the general situation in 
which the distribution P of 5n is merely known to belong to a 
set of distributions V. Specifically, we have that V correspond 
to {N{0, 1)} in the Gaussian case 1) and to {ti, ^0, \/('^^^)) } 
in the Student’s t case 2) with ly degrees of freedom. In the 
former two cases, the set 7^ is a singleton as the distribution 
of Sn is known. In the most general case 5), the set P 
consists of all distributions of zero mean and unit variance. 
The set P is additionally required to contain only unimodal 
or unimodal symmetric distributions in the unimodal case 4) 
and the unimodal symmetric case 3), as both notions are scale 
invariant. That is, if 5 is unimodal or symmetric unimodal than 
a6 + /3 is unimodal or symmetric unimodal as well for any a 
and (3 real numbers. 

Given a set of distributions P representing the distributional 
assumptions made, we will first show that the chance con¬ 
straint (fTSl l is equivalent to the deterministic constraint (fT4l) 
for /•p(fc) := infpep P [(5„ < k]. 

We can trivially rewrite constraint (l20l) equivalently as the 
constraint P[(5„ < (c — /r(PG))/ o'(P’g)] > 1 — e for all P € "P. 
Using the definition of f-p we can demand alternatively that 
fv{{c — fi^Pg)) /(j{Pg)) > 1—e. Einally as the function fp is 
increasing it has a well defined generalized inverse fp^{X) = 
inf {k I fp{k) > A}. Prom the definition of / it follows that 
we must have that (c — ^,{Pg)) /<^{Pg) > /^^(l ~ e)- After 
reordering the terms we obtain 


a{PG) + biPG)iiR < c - f^\l - e)\\biPG)X:](^h. 


Einally in the remainder of this section, we show for 
the distributional assumptions discussed above we obtain the 
results mentioned in Table U 

1-2) Por non-atomic distributions P it follows by continuity 
that /{pj. = P[(5„ < fc] = P[5n < k] explaining the first 
two entries in Table |I] 


















3) Let the set V corresponds to the set of all symmetric uni- 
modal distributions S with zero mean and unit variance 
then we can leverage the classical Gauss bound 1^ . 
Il24l . Indeed, in this situation we have that the mode of 
the distributions is known as it coincides with the mean 
because of symmetry. For any fc > 0, we have again by 
symmetry that the equality P[(5„ > k] = iP[|(5„| > k] 
holds. As the mode of P is known, the classical Gauss 
bound ED, in can be used to establish 


fu{k) 


1 — sup P [(5„ > k] 
pgs 


r 1 - isuppgy P[|^„| > k] 

I 0 


9fc^-2 

9fe^ 



if fc > 
fc > 0 
otherwise 


if fc > 0 
otherwise 


( 21 ) 


4) The unimodal case, in which V = U consists of zero 
mean and unit variance unimodal measures, can be 
dealt with using the one-sided Vysochanskij-Petunin 
inequality m . i.e. 


fu{k) = 


1 — sup P [5n > fc] 

PG(7 


= 1 - 


4 

9(T+W) 


if fc > 


1 - 

4 

3 


1-It^ iik>o 

1 

otherwise 

4 

9(l + fc^) 

if fc > 

fc" 

l + fc2 

fc > 0 

0 

otherwise 


( 22 ) 


5) Lastly, when the set V corresponds to the set C con¬ 
taining all distributions of zero mean and variance, the 
classical Cantelli inequality m establishes that 


/c(fc) = 1 - sup P[(5„ > fc] 

pgc 


= 1 - 


1 

l-|-fc2 

1 


fc" 

14-fc" 

0 


if fc > 0 
otherwise 

if fc > 0 
otherwise 


(23) 


which after taking the inverse gives us the corresponding 
result in Table U 
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